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ABSTRACT 

The cosmic star formation rate (CSFR), namely the star formation rate in a unitary 
comoving volume of the Universe, is an important clue to investigate the history of 
the assembly and evolution of galaxies. Here, we develop a method to study the CSFR 
from a purely theoretical point of view. Starting from detailed models of chemical 
evolution, which best fit the properties of local galaxies, we obtain the histories of 
star formation of galaxies of different morphological types (ellipticals, spirals, irregu- 
lars). These histories are then used to determine the luminosity functions of the same 
galaxies by means of a spectro-photometric code. We obtain the CSFR under different 
hypotheses about galaxy formation scenarios. First, we study the hypothesis of a pure 
luminosity evolution scenario, in which all galaxies are supposed to form at the same 
rcdshift and then evolve only in luminosity without any merging or interaction. Then 
we consider scenarios in which the number density or the slope of the luminosity func- 
tions are assumed to vary with redshift. After comparison of our results with the data 
available in literature, we conclude that a pure luminosity evolution does not provide a 
good fit to the data, especially at very high redshift, although many uncertainties are 
still present in the data because of the unknown dust corrections and assumed initial 
mass function. On the other hand, a variation in the number density of ellipticals and 
spirals as a function of redshift can provide a better fit to the observed CSFR. We 
critically discuss the possible scenarios for galaxy formation derived from this finding. 
We also explore cases of variable slope of the luminosity functions with redshift as 
well as variations of number density and slope at the same time. We cannot find any 
of those cases which can fit the data as well as the solely number density variation. 
Finally, we compute the evolution of the average cosmic metallicity in galaxies with 
rcdshift and show that in the pure luminosity evolution case the ellipticals dominate 
the metal production in the Universe, whereas in the case of number density evolution 
are the spirals the main producers of metals. 

Key words: galaxies: evolution - galaxies: fundamental parameters - galaxies: pho- 
tometry - galaxies: luminosity function - ISM: general - ISM: evolution. 



1 INTRODUCTION 

The rate at which galaxies have formed stars throughout the 
whole cosmic time is a fundamental clue to investigate the 
history of the assembling and evolution of structures in the 
Universe. The cosmic star formation rate (CSFR), defined 
as the comoving space density of the global star formation 
rate in a unitary volume of the Universe, is not a directly 
measurable quantity. It can be inferred only from the mea- 
surement of the luminosity density in different wavebands, 
which are then transformed into star formation rate by suit- 
able calibrations. Starting from its first determination made 
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bv lLillv et all (|l996l l and thanks to deeper and deeper sur- 
veys it has been possible to extend the d etermination of the 
cosmic star formation rate up to z ~ 10 jMadau et al.l 19981: 
Hopkins! 120041: IBouwens et all l200Sl: 120081: IKistler et al l 



20091 : IBouwens et alj|2011al : Tlshida et al.ll201lf ). Constrain 



ing the CSFR at high redshift would be of paramount impor- 
tance to put constraints in the galaxy formation scenario. In 
particular, to decide between the Monolithic Collapse (MC) 
and the Hierarchical Clustering (HC), because the CSFR de- 
pends on both the star formation rate (SFR) in galaxies and 
on the evolution of the galaxy luminosity function (LF). In 
the MC scenario, spheroids and bulges formed at high red- 
shift (e.g. z > 2 — 3) as the result of a violent burst of star 
formation, following the collapse of a gas cloud. After the 
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development of a galactic wind quenching st ar formation , 
galaxies evolved passively un til present ti me (lLarsonlll974l ; 
IMatteucci fc Tornamblll987l ). Moreover. [Matteuccil l| 19941 ) 
introduced in this scenario the assumption that more mas- 
sive spheroids had higher star formation efficiencies than less 
massive ones, and stopped to form stars earlier. It is common 
now to refer to this behaviour as "downsizing in star forma- 
tion" . These assumptions allow us to reproduce the majority 
of chemical and photometric properties of local ellipticals. 
On the other hand, in the HC scenario for galaxy formation, 
baryon assembly basically mirrors the hierarchical build up 
of dark matter structures. Therefore, spheroids and bulges 
formed as a series of subsequential mergers among gas-rich 
galaxies or with galaxies that have already stopped their star 
formation. In this scenario, galaxies form stars at lower rates 
than in the MC scenario, with more massive s pheroids reach- 
ing their final mass at later times (z < 1.5) (White fc Reed 
1 19781 ; iBaugh et al.lll998l ; ICole et al.ll2000l ; iMenci et al.ll2004h . 
In this paper, we calculate the evolution of the CSFR and 
the mean metallicity of the interstellar medium (ISM) in 
galaxies. This is done by means of detailed chemical evolu- 
tion models for galaxies of different morphological types, i.e., 
ellipticals, spirals, and irregulars, which successfully repro- 
duce the local properties of such galaxies, and from which we 
derive their star formation histories (SFH). Star formation 
histories obtained with chemical evolution models are then 
used to determine the photometric evolution of the galaxies 
(spectra and magnitudes), through a spectro-photometric 
code. This allows us to compute the LF, that gives the num- 
ber of galaxies in a unitary volume of the Universe in a lumi- 
nosity bin. In t his work, we normaliz e our LFs to the local 
LFs derived by iMarzke et al.l (|l998T l. which are based on 
the "Second Southern Sky Redshift Survey" (SSRS2) using 
data from 5404 galaxies at z 0.05. These LFs are deter- 
mined in the B band for galaxies of different morphological 
type. Then we compute the evolution in luminosity of the 
galaxy at the break of the LF for each galactic type. This is 
done under different scenarios: first, we study the hypoth- 
esis of a pure luminosity evolution scenario (PLE). In this 
case, galaxies are supposed to form all at the same red- 
shift and then evolve only in luminosity without any merg- 
ing or interaction. In this case, we follow the approach of 
ICalura fc Matteuccil |2003l ). Then we define two new scenar- 
ios: one in which the number density of the LF is assumed 
to vary with redshift, and one in which the slope of the LF 
is assumed to vary with redshift. All the results have been 
compared with the available data found in literature. Fi- 
nally, we calculate the mean metallicity of the ISM for the 
PLE and the number density evolution (NDE) scenario. The 
paper is organized as follows: in § 2 we describe the chemical 
and photometric models, in § 3 we present our results and 
in § 4 we draw some conclusions. 

Throughout this paper the cosmology adopted will fol- 
low the KCDM paradigm with ilo = 0.3, Q.\ = 0.7 and 
h = 0.65. 



2 CHEMICAL AND PHOTOMETRIC 
EVOLUTION OF GALAXIES 

Chemical evolution models allow us to compute the evolu- 
tion in time and space of several quantities such as the star 



formation, the production rate of chemical elements and the 
chemical abundances in the stars and in the gas, starting 
from the matter reprocessed by the stars and restored into 
the ISM through stellar winds and supernova explosions. 
Starting from observational constraints such as the present 
day chemical abundances, it is possible to perform a detailed 
backward evolution of galaxies of different morphological 
type. We consider three galactic types namely irregulars, 
spiral disks and spheroids which include both bulges and el- 
liptical galaxies. This is motivated by the fact that ellipticals 
and bulges are likely to have a common origin, i.e. both are 
likely to have for med their stars on ve ry short timescales and 
a long time ago ( Jablonka ct al. 1996)) . Detailed descriptions 
of the adopted mo dels can be found in Matteucci fc Francois! 
(1 19891) for spirals. IPipino fc Matteuccil (|2004h for ellipticals 
and Bradamante et al. ( 19981 ) for irregulars. 

In our picture, spheroids form as the result of a sudden 
collapse of a gas cloud of primordial chemical composition. 
The gas is supposed to be accreted over a finite timescale, 
so in this respect we do not assume pure MC models. A 
galaxy can be modelled either as a one-zone object or a 
multi-zone object in which the galaxy is divided into several 
concentric shells independent of each-other. The spheroids 
are modelled as a one zone objects and they are supposed 
to suffer a strong initial burst of star formation, that stops 
as soon as the galactic wind develops; in this sense their 
evolution crucially depends on the time at which the galactic 
wind occurs. From that moment on the galaxy is supposed 
to evolve passively. 

The chemical evoluti on of spiral disks is studied b y us- 
ing the one-infall model of lMatteucci fc Francois! (| 19891 ). de- 
veloped to reproduce the observational constraints of the 
Milky Way disk. This model belongs to the class of the so- 
called "infall models", where the disk is thought as gradually 
forming by infall of primordial gas. Instantaneous mixing of 
gas is assumed, but the instantaneous recycling approxima- 
tion is relaxed. The galactic disk is supposed to have formed 
on timescales increasing with the galactocentric distance, 
being larger at larger radii, in agre ement with th e dissipa- 
tive models of galaxy formation of iLarsonl (1976). The gas 
accumulates faster in the inner than in the outer region, 
according to the so-called "inside-out" scenario. The pro- 
cess of disk formation is much longer than the halo and the 
bulge formation (they should form on timescales no longer 
than 1 Gyr) with typical timescales varying from ~ 2 Gyr 
in the inner disk, ~ 7 Gyr in the solar region and up to 
15 — 20 Gyr in the outer disk. This mechanism is impor- 
tant to reproduce the observed abundance grad ients (see 
IMatteucci fc Francois|[l989l , IChiappini et al.ll200lh . The spi- 
ral disk is approximated by several independent rings, 2 kpc 
wide, without exchange of matter between them. In this 
work the disk has been set to extend from 2 to 16 kpc. 

Irregular galaxies are modelled as one zone objects and 
are assumed to assemble all their gas by means of a contin- 
uous infall of primordial gas and to form stars at a lower 
rate than the other morphological types. Their stellar pop- 
ulations appear to be mostly young, their metallicity is low 
and their gas content is large. All these features indicate 
that these galaxies are poorly evolved objects. Furthermore, 
irregular galaxies are allowed to develop moderate galac- 
tic winds which do not halt the star formation as in the 
spheroids l|Bradamante et alii 19981 ). 
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2.1 The chemical evolution models 

The equation that describes the evolution of the i-th element 
has the following form: 
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Where Gi(t) is the normalized fractional mass of gas 
within a galaxy in the form of the element i: 

M g (t)Xj{t) 



Gi{t) = 



Mtot 



(2) 



where Mtot and M g (t) are the total galaxy mass and the 
mass of gas at time t, respectively. 

Gi(t) 



Xi(t) = 



G(t) 



(3) 



represents the abundance in mass of the element i, the sum- 
mation over all elements in the gas mixture being equal to 
unity. Thus, the quantity: 



G(t) 



Mg(t) 

Mtot ' 



(4) 



is the total fractional mass of gas present in the galaxy at 
time t. In the case of spiral disks M g (t) and Mtot are replaced 
by the corresponding surface densities a(r, t) and a(r, to) 

The first term on the right hand side of eq. [TJ including 
the star formation rate tp(t), represents the rate at which 
each element disappears from the ISM owing to the star 
formation. 

The second term is the rate at which each element 
is restored into the ISM by single stars with masses in 
the range [Ml, Mb™.], where Ml is the minimum stellar 
mass contributing, at a given time t, to chemical enrich- 
ment (~ O.8M0) and MBm is the minimum total binary 
mass of systems which can give rise to Type la SNe (3 Mq, 
iMatteucci fc Greggidll986T ). The quantity Q m i(t — r m ) indi- 
cates the fraction of mass restored by the stars in form of an 
element i that can be produced or destroye d in stars or both. 
This is the so-called "production matrix" l|Talbot fc Arnettl 

[lizi). 

The third term represents the enrichment due to bi- 
naries which become Type la SNe, i.e. all the binary sys- 
tems with total mass in the range between Mg m and Mbm, 
with Mbm being the maximum total mass of a binary sys- 
tem which can give rise to Type la SNe and is equal to 
16 Mq (8 + 8 Mq). The model adopted for the progeni- 
tors of the Type la SNe is the single degenerate one as in 



IMatteucci fc Greggiol (jl986T ). The parameter A represents 
the fraction of binary stars giving rise to Type la SNe, and 
although its real value is unknown it is fixed in order to re- 
produce the observed present time SN la rate. This param- 
eter obviously depends on the adopted IMF and, in general, 
values between 0.05 and 0.09 are adopted. This ensures to 
reproduce the actual SN la rate both in the Milky Way and 
in other galaxies. It is worth noting that, in this term, the 
quantities ip and Q m i refer to the time (t — T m 2), where r m 2 
indicates the lifetime of the secondary star of the binary sys- 
tem, which regulates the explosion timescale of the system. 
The coefficient fi = M2/MB, is the ratio between the mass 
of the secondary component and the total mass of the binary 
system, while /(/i) is the distribution function of this ratio. 
It is calculated on a statistical basis indicating that values 
close to 0.5 are preferred. Its analytical expression can be 
written as: 
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with 7 = 2 as a parameter (|Greggio fc Renzinil [l983h and 
Umin is the minimum mass fraction contributing to the SNIa 
rate at the time t, and is given by: 
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The fourth term of eq. Q] represents the enrichment due to 
stars in the range [Msm, Mbm] which are single or, if bina- 
ries, do not produce a SNIa event. All the stars in this range 
with mass larger than 8Mq are assumed to explode as core 
collapse supernovae. 

The fifth term represents the contribution to the chem- 
ical enrichment coming from stars more massive than Mbm 
(supposed to explode like Type II and l b / c SNe). In all the 
models the upper mass limit contributing to chemical en- 
richment is assumed to be 100 Mq. 

Finally, the last two terms are the rate of accretion of 
matter with primordial abundances Xa, and the outflow 
rate for the element i, respectively. 



2.1.1 The Star Formation Rate 

The main feature characterizing a particular morphological 
galactic type is represented by the prescription adopted for 
its star formation history. In the case of elliptical galaxies 
the SFR ijj has the simple form given by: 



yj) = uG{t), 



where v is the star formation efficiency, namely the inverse 
of the typical time scale for star formation, and assumed to 
be 15 Gyr -1 . 

In the case of irregular galaxies, a continuous star for- 
mation rate is assumed as in eq. [7J but characterized by 
a lower efficiency than the one adopted for ellipticals (0.1 
Gyr -1 in this case). 

In the case of spiral galaxies the fo rmu lation is the fol- 
lowin g (see IMatteucci fc Francoidll989l and lChiappini et"al1 
Il997l) : 



ip(r,t) = v 
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Morphological 


Mtot 




V 




Type 


(Mq) 


(Gyr) 


(Gyr- 1 ) 


(Kpc) 


Ellipticals 


10 11 


0.3 


15 


3.0 


Spirals 


4xl0 10 


7.0 


1.0 


3.5 


Irregulars 


10 10 


0.5 


0.1 


1.0 



Table 1. Model parameters for the galaxies reproduced in this 
work. Mtot is the baryonic mass, R e ff is the effective radius, r in j 
is the infall timescale and v is the star formation efficiency. For 
spirals we indicate the average infall time scale for the disk. 



where v is set in this case equal to 1 Gyr -1 , o~(r,t) is the 
total surface mass density at a given radius and at a given 
time, o~(tq, i) is the total surface mass density at the position 
of the Sun (set equal to 8.0 Kpc). Finally a(r,ta) is the 
surface mass density at present time . The exponent k is 
assumed equal to 1.4, as suggested by iKennicuttl l|l998T ). 

All these models are tuned to reproduce at best the 
main observed properties of local galaxies of different mor- 
phological type. In Table[T]we show the main characteristics 
of the galaxies of different morphological type that we use 
in this work. It is important to stress that in our models, we 
consider only one galaxy per morphological type and take 
it as representative of the whole population. For the spiral 
galaxy we take the average SF and metallicity of the disk. 

In Figure Q] we show the histories of star formation of 
the g alaxies of d ifferen t morphological type described above, 
for a ISalpeteij (jl955T ) IMF. It is worth noting that for a 
spiral like the Milky Way an IMF stee per than t he on e used 
is usually adopted. Here we assume ISalpete 3 ll 19551 ) IMF 
for all galaxies because the CSFR derived observationally is 
obtained by means of this IMF. As it is clear from the Figure 
a typical elliptical (10 11 Mq of luminous mass) experiences 
a high burst of star formation lasting for about 0.6 Gyr with 
a maximum peak of more than ~ 100 M^yr' 1 . After such a 
period, the star formation ceases abruptly owing to the onset 
of the galactic wind. The spiral disk is characterized by a 
continuous SFR with a large peak around 1 ~ 2 Gyr of less 
than ~ 10 M0i/r _1 and a present time value of ~ 2 Moj/r -1 
(for the Milky Way). Finally, the irregular galaxy forms stars 
at a rate smaller than the previous morphological types, 
with an increasing SFR that reaches a maximum of less then 
1 M0j/r _1 , with a smaller present value, due to the onset of 
the galactic wind. 



2.2 The photometric model 

To study the evolution of the photometric properties 
of galaxies, we u se the spectro-photometric code GRASIL 
(|Silva et al.l I1998T ). Its starting point is the output of a 
chemical evolution code that follows the time evolution of 
SFR, metallicity (Z) and gas fraction. The building blocks 
of galaxy photometric models are the library of isochrones 
and the stellar atmosph eric models. GRASIL is based on the 
Padova stellar models (jBertelli et al.|[l994l ). with a major 
difference consisting in the computation of the effects of 
dusty envelops arou nd AGB stars , and on the stellar at- 
mospheric models by iKurucdj 19921 ). The single stellar pop- 
ulations (SSPs) span a range from 1 Myr to 18 Gyr for what 
concerns ages, and Z = 0.004, 0.008, 0.02, 0.05, 0.1 for the 
metallicities (keeping the relative proportion of the metals 




Time [Gyr] 



Figure 1. Star formation rates as a function of time, for a typical 
elliptical of 10 11 Mq (green solid line), spiral of ~ 4 X 10 10 Mq 
(red dotted line) and irregular of 10 10 Mq (blue dashed line) 
galaxy. 



equal to the solar partition). At this point, the spectral syn- 
thesis is performed summing up the spectra of each stellar 
generation provided by the SSP of the appropriate age and 
metallicity, weighted by the SFR (ip) at the time of the star 
birth. The flux at a given wavelength is defined as: 



F x (t G )= ip(t)xSSP x (t G -t,Z(t))dt, 



(9) 



where Fx{ta) is the integrated monochromatic flux at the 
present time Ig and SSP\(tc — t,Z(t)) is the integrated 
monochromatic flux of an SSP at age t and metallicity Z. If 
the effects of dust are neglected (as in our case) the spectral 
energy distribution (SED) is simply given summing up the 
spectra of all stars. 



2.3 The luminosity evolution of galaxies 

To derive the CSFR, it is of fundamental importance the 
knowledge of the relative distribution of galaxies, of any 
morphological type, as a function of redshift. This can be 
done through the LF, which gives the distribution of galax- 
ies per unit volume in the lumin osity interval \L , L + dL]. 
The LF is well reproduced by the lSchechterl (|l976l ) function: 




dL 

L*' 



(10) 



where <f>* is the number of galaxies per unit volume; L* is 
the characteristic luminosity which separates bright sources 
from faint sources and a is the slope of the luminosity 
function. Different determinations of the local LF exist (see 
Ide Lapparentl l2003 for a review) and they show several dif- 
ferences, in the sense that the parameters of the Schechter 
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Morphological 


M* 


a 




Type 








Ellipticals 


-19.37 


-1.00 


4.4 


Spirals 


-19.43 


-1.11 


8.0 


Irregulars 


-19.78 


-1.81 


0.2 



ca n calculate the CSFR acco rding to the method adopted 
bv lCalura fe Matteuccil l|2003r ): 



Table 2. Parameters of the LF as defined in Marzkc ct al. (1998). 
Here M* refers to the magnitude of the galaxy at the break, <f>* 
refers to the normalization of the LF (expressed in Mpc~ 3 ) and 
a is the slope of the LF. 



function are not fully constrained by observations. The lu- 
minosity function can be measured in several bands, obvi- 
ously depending on the redshift regime under investigation. 
In this work, we make use of th e determination of t he lo- 
cal B-band luminosity function of lMarzke et al] l|199a ). that 
has the characteristic of considering separately the different 
galactic morphological types. For each morphological type 
the parameters of the LF are summarized in table [21 

Once we have the LF, we can compute the luminosity 
density (LD), which is what is really observed. The lumi- 
nosity density, in a given band, is defined as the integrated 
light radiated per unit volume from the entire galaxy popu- 
lation. It stems from the integral over all luminosities of the 
observed luminosity function: 



PL = 



dL. 



(11) 



At z — 0, the LDs for the single galaxy types are given by 
the above integral of the present time luminosity functions of 
iMarzke et al] (1998). At redshifts other than zero, for each 
morphological type we consider the luminosity evolution ob- 
tained with the spectrophotometric code. 



3 RESULTS 

As already pointed out, the CSFR is not a directly measur- 
able quantity. It can be inferred only from the measurement 
of the luminosity density in different wavebands and then 
converted into SF by means of suitable calibrat ions. Its first 
determination was made by iLillv et al] l|l99rj) (at 2800 A, 
4400 A, and 1 u rn), followed in the subseq uent years by 
iHogg et al] <| 19981 ) (using the [Oil] line) and IWilson et all 
( 20021 ) (in the UV). These studies constrained the CSFR 
from z = to z = 1, establishing that it is steadily increas- 
ing in that range. The CSFR seems to peak around z ~ 2 — 3 
and then decrease at higher reds hifts, where the r e is n ot a 
well defined trend. For example iBouwens et al. 1 20081) on 
the basis of optical and UV data and lHopki~ ( 20041 ) on 
the basis of a compilation of data taken in different bands, 
found that the CSFR from roug hly z ~ 3 to z ~ 9 is steadily 
decreasing. On the other hand. iKistler et all l|2009l ). on the 
basis of the gamma ray burst (GRB) data from Swift found 
a CSFR that is roughly constant from z ~ 4 to z ~ 8. 



3.1 The cosmic star formation rate 

By means of the star formation histories of galaxies of differ- 
ent morphological type and their photometric evolution we 




(z)ipi(z), 



(12) 



where p*(z) is the SFR density, ps% is the B-Band luminosity 

density, ^irj ^ s l ne B-band mass-to-light ratio and ipi 

(expressed in yr~ x ) is the star formation rate per unit mass 
of the i-th morphological type. 

In this work, we make use of the B-band luminosity, 
because it is a good star formation tracer since hot massive 
stars emit in this band. 

In principle, the CSFR, could be computed simply by 
summing the SFR of each galaxy type multiplied by the cor- 
responding number density. However, this simple approach 
would not allow us to explore in detail the case in which the 
slope of the LF varies with redshift (see section 13 . 1 . 3 p . We 
tested that the two approaches lead to the same result. 



3.1.1 The Pure Luminosity Evolution Scenario 

Here we describe the determination of the CSFR in the 
framework of a PLE scenario. With pure luminosity evo- 
lution, we mean that galaxies of all morphological types are 
supposed to form at high redshift and then to evolve only 
in luminosity and not in number. We also assume that the 
slope of the LF is constant with redshift. In this context, 
galaxies are assumed to be isolated, i.e. the effects of merg- 
ers or interactions are irrelevant at all redshifts, and it is 
assumed that they start forming stars all at the same time. 

In this simple approach, the differences between the var- 
ious morphological types should be ascribed only to their 
different star formation histories. Early-type galaxies (i.e. 
ellipticals, bulges and SO) are dominated by old stellar pop- 
ulations so they must have formed their stars with high ef- 
ficiency and in a short period, while late-type galaxies (i.e. 
spirals and irregulars) should have formed stars for the whole 
Hubble time, as seen in Fig. [T] 

In Figure [2] we can see the predicted evolution of the 
luminosity density for ellipticals, spirals and irregulars, in 
the U (centered at 3650 A), B (centered at 4450 A), I (cen- 
tered at 8060 A) and K (centered at 21900 A) band. These 
evolutions are computed by adopting the SFRs of Fig. [T] 
and constant slope and number density in the luminosity 
function. 

From the Figure it is clear that, at early times in ev- 
ery band, the luminosity density is dominated by the light 
emitted by elliptical galaxies. This is obviously due to the 
fact that in the monolithic scenario they all suffer a strong 
initial burst of star formation which lasts for < 1 Gyr and 
then evolve passively. At this time a galactic wind develops 
and star formation stops, so their luminosity in the U and B 
bands, where young, newborn stars emit, begins rapidly to 
decrease. On the other hand, spirals form stars continuously 
and the decrease in luminosity in the U and B band should 
be ascribed mostly to the consumption of their gas reservoir. 
Finally, also irregulars form stars for the entire cosmic time 
but their luminosity density is lower than that of the other 
morphological types due to their low star formation. Their 
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time [Gyr] time [Gyr] 



Figure 2. Prediction of the luminosity density evolution for 
galaxies of different types in U (top left), B (top right), I (bottom 
left) and K (bottom right) band respectively, in the case of pure 
luminosity evolution. Here the x-axis represents the elapsed time 
since z=10. Solid green line refers to ellipticals, red dotted line to 
spirals and long dashed blue line to irregulars. 



luminosity density is also decreased by the onset of galactic 
winds at late times. 

The I and K bands are dominated by the light emitted 
by less massive, long living stars. Also in this case, ellipti- 
cals are the main sources at early times, while spirals become 
dominant for cosmic time > 2 Gyr. However, the difference 
in the values reached by ellipticals and spirals is less pro- 
nounced than in U and B bands. The luminosity density 
of spirals, after an initial phase of constant increase lasting 
~ 3 Gyr is observed to decrease slowly in the I band and to 
be quite constant in the K band. Irregulars, also in this case 
have values lower than ellipticals and spirals which tend to 
decrease very slowly through the whole cosmic time. 

The PLE model belongs t o the category of the so -called 
"backward evolution" models l|Tinslev fc Danlvlll976l ); these 
models start from a well known and well constrained de- 
scription of the local Universe to reconstruct the evolution 
of galaxies back to their formation. 

In our model, all the galaxies are supposed to form at 
z = 10. This choice is motivated by the fact that more 
and more objects are re vealed at z > 5 (Mobashc r et al.l 
120051 IVanzella e t al. 2008), with evidences of a z ~ 10 detec- 
tion (jBouwens et alj|2011bl ). The predicted CSFR is shown 
in Figure [3] for compariso n we plot in the same Figur e also 
the CSFR as obtained by ICalura fc Matteuccil (|2003l ) who 
adopted very similar assumptions and star formation histo- 
ries. 

As it can be seen, both models show some discrepancies 
when compared with the observational data. In our case, the 
discrepancy begins at z ~ 0.5 and it is particularly evident 
in the redshift interval running from ~ 0.5 to ~ 5. Here 



Figure 3. Cosmic star formation history in the case of pure lu- 
minosity evolution. Blue solid line, model described in th is work; 
green dashed li ne model from | Calura fc Mattcucci ( 2 0031). Points 
are taken from: | Hop_kirj jgOOj) red open circles, iLv et alj (120071 1 
bro wn sq uares, iTakahashi et all d2007h purple open squares, 
cyan open h e xagon s, ISevmour et al l [120081) yellow 
pentagons, iKistler et alj J2009I) ora n ge cros ses, ILv et al.l boill ) 
magenta triangles, Bouwens et al.l d201 laT) green open stars, 
ICucciati et"afl ll201lh light green crosses. 

we underpredict the observed values by a factor of ~ 3. 
This is the redshift interval in which only spirals and irreg- 
ulars contribute to the cosmic star formation history. Also 
the Calura & Matteucci CSFR suffers of the same problem 
when compared with data. On the other side, at high red- 
shift, our model seems to overpredict the data, as well as 
that of Calura & Matteucci. The high observed peak is to 
ascribe to ellipticals which, in order to reproduce the local 
observational features, (like the increase of the [Mg/Fe] ra- 
tio with galactic mass) have to form stars in a short time 
interval (lasting less than 1 Gyr, lMatteucci|[l994l ). This im- 
plies that the SFR of these objects must lie in the range 
~ 100 - 1000 Mcyr' 1 . At z ~ 6 ellipticals stop to form 
stars and the CSF R decreases abruptly . The d iscrepancies 
with the model of ICalura fc Matteuccil (|2003l ) can be ex- 
plained through the differences in the parameters adopted 
in the chemical evolution mod els and in the different pho - 
tometric code used. In fact, in | Calura fc Matteuccil (2003), 
the photometric code of ljimenez et al.l (|l998l ) was adopted. 
Moreover, the chemical mo del for the ellipticals adopted by 
ICalura fc Matteuccil (|2003D is a closed box one, whereas our 
model includes gas infall. 

3.1.2 The Number Density Evolution Scenario (NDE) 

In this section, we describe a different approach for the calcu- 
lation of the CSFR. In what follows we introduce a scenario 
in which we let the number density of galaxies to evolve with 
redshift while the SFHs of galaxies are the same of Fig. [T] 
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The main consequence of relaxing the hypothesis of 
the constancy of the number density of galaxies is that 
we consider galaxy interactions. It has been claimed that 
giant ellipticals could form, as the result of the occur- 
rence of 1 to 3 major dry -mergers during galactic lifetime 
ijPipino fc Matteuccil 120081 and reference therein) and this 
view is in agreement with the observations. Moreover it 
has been observed that the fraction of dry-mergers i ncreases 
with redshift from z ~ 1.2 to z = (|Lin et aljfcoQgh . 

To include the number density variation we modify the 
expression of the luminosity function as defined in eq. 1101 
letting the parameter <f)* (i.e. the number density of galax- 
ies) to vary with the redshift. The new formulation of the 
luminosity function is: 

*(£)§ = Mi + z { - L/L ' ] ^ (is) 

where 0o = 4>*{ z = 0), i.e. the present time number density 
of galaxies. This model is constrained in orde r to reproduce 
the nu mber densities observed at z = by iMarzke et al.l 
l|l998h . 

We let the parameter ft to run from —6 to +6 in steps 
of 0.2. The choice of this interval has been made in order to 
explore a set of values as wide as possible without giving us 
meaningless results and, at the same time, letting us to test 
cases of extreme number density variation. For each value 
of ft, we calculate the luminosity density using the modified 
version of the Schechter function. Then, using the star for- 
mation histories depicted in section [2. 1.1 1 we determine the 
star formation rate density, using the equation 1121 

We perform a careful check of the contribution to the 
total CSFR from the different morphological types, and on 
its basis we decide to let only spirals and ellipticals to evolve. 
Therefore, we consider the number density of irregulars con- 
stant through the cosmic time. The reason resides in their 
contribution to the total CSFR, that is, at any redshift, two 
orders of magnitude lower than the one of ellipticals and 
spirals and therefore negligible. Our best values of the ft 
parameter for each morphological type are: 

• Ellipticals: ft = -0.8; 

• spirals: ft = 1; 

• irregulars: ft = 0. 

This means that the number density of ellipticals is sup- 
posed to increase from z — 10 to z — (with ~ 15% of ellip- 
ticals already in place at z = 10). The opposite behaviour 
is predicted for spirals, which are supposed to decrease with 
decreasing redshift. Under these assumptions, there is a pos- 
sible interpretation in terms of galaxy formation. For exam- 
ple, the predicted number density evolution could mean that 
a large fraction of elliptical galax ies can form thanks to dry 
mer gers of spirals, as pro] josed bv lToomre fc Toomrel (|l972l ) 
and iNaab fc Burkertl (|2003m . Since we are considering only 
dry-mergers it means that the merging process quenches the 
star formation in the progenitors and that the residual gas 
is lost. From our model we can see that the number density 
of spirals decreases of roughly 60% from z ~ 2 to z ~ 0, 
and of 50% from z ~ 1 to z = 0. This is in good agreement 
with iBoissier et all (|2010h who found a decrease of 50% of 
Milky Way siblings from z ~ 1 to z = 0, indicating this 
percentage as a lower limit. In Figure [4] we show our best 



model for the CSFR in the case of number density evolution. 
From the plot we can see that the peak observed in the PLE 
scenario completely disappears. The model seems to better 
reproduce the high redshift trend if the error bars are taken 
into account. Also in the redshift interval between and 
~ 5 the agreement is improved although the curve is still 
slightly lower than the data, but this is a problem common 
to all the other models of CSFR, as we will see later. Also 
in this case we underpredict the decrease of the CSFR be- 
tween z ~ 2 and z — 0. This, in principle, can be due to the 
lack in our models of strong starburst galaxies at low and 
intermediat e redshift (see e.g . lElbaz et ai]|201ll ). even if as 
claimed by iRodighiero et alT l|201ll ) starbursts seem to ac- 
count only for the 10% of the SFR density at z ~ 2. On the 
other hand, our predicted slope of the C SFR in this redshift 
range good as well as for the model of Ivan de Voort et all 
(|20 1 IT) . These authors using cosmological smoothed particle 
hydrodynamics simulations, explain this drop with a decline 
in the cold-accretion rate density on to haloes and with AGN 
feedback that decreases the contribution to the CSFR of the 
gas accreted in the hot mode. 

In Figure [5] we show the predicted evolution of the lu- 
minosity density, for our best model, in the case of num- 
ber density evolution. Comparing these predictions with the 
ones in Figure [2] it is possible to see that in this case spirals 
dominate the luminosity density in every band and at any 
epoch, at variance with the PLE case. 

It must be said that the variation of the number density 
of galaxies of different morphological types is still matter of 
debate, since there is, in ge neral, no concordance among 
various authors. For example, iTotani fc Yoshiil (|l998h found 
that, parametrizing the comoving number density of E/S0 
as proportional to (1 + z) 7 , the favoured value for 7 is — 0.8± 
1.7, in perfect agreement with our value. This analysis is 
valid in the redshift range 0.3 < z < 0.8, and the large 
error is due to the limite d redshift range of their sample. 
Ide Lapparent et all (|2004h , using data from the ESS (ESO- 
Sculptor Survey), found a number density evolution of late- 
type galaxies following the same parame trization, w it h 7 = 
2 ± 1. On the other hand, for example iLillv et all (Il998l ) 
found at z < 1 an evolution of large disk structures (with 
a bulge-over-total ratio B/T > 0.5) following a (1 + z) 1 
relation, with 7 = ±0.5, lower than the one found in our 
model. 

Even more challenging is to asses the evolution at high 
redshift, since very deep surveys are required and the cor- 
respondence between actual galaxies and their high z coun- 
terparts is not firmly established. This causes a loss of infor- 
mation about the morphology of the various objects giving 
controversial results. Moreov er, the uncerta in dust correc- 
tions and the adopted IMF (|Salpeterl fl955'l used to derive 
the CSFR, could underestimate the actual CSFR. There- 
fore, we do not completely exclude that the CSFR could be 
constant or increasing at high z. 

3.1.3 The Variation of a 

In this section we describe a modification of the Schechter 
function in which we consider the evolution with redshift of 
the slope of the LF for galaxies of different morphological 
type. We start again from the LF as defined in sec. 12.31 and 
assume that a — aa(l + z) 13 , where ao = a(z — 0) is the 




Figure 4. The cosmic star formation history in the number den- 
sity evolution scenario. Blue line refers to model results. Data 
points are the same described in Figure [3] 



Figure 6. The CSFH in the case of variation of the LF slope a. 
Here we show the cases corresponding to /3 = 0.4 (top panel) and 
ft = —0.2 (bottom panel). Blue lines indicate the results of the 
model. Data points are the same described in Figure [3] Note the 
different y-axis scale adopted for clarity. 
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Figure 5. Prediction of the luminosity density evolution for 
galaxies of different types in U (top left), B (top right), I (bottom 
left) and K (bottom right) band respectively in the case of num- 
ber density evolution (best model). Here the x-axis represents the 
elapsed time since z=10. Solid green line refers to ellipticals, red 
dotted line to spirals and long dashed blue line to irregulars. 



slope of the LF at present time and (3 is a free parameter. 
Therefore, the LF assumes the form: 



*(L) 



clL 



dL 



(14) 



This means that at any redshift we change the slope of 
the LF. As in the case of the number density evolution, we 
let /3 to run from —6 to +6 in steps of 0.2, and in no case 
we have found an improvement in fitting the data relative to 
the pure luminosity or number density variation scenarios. 
From the observational point of view, there are some sug- 
gestions of variations of the slope of the LF but referred to 
the whole population of galaxies, w ithout differentiat i ng for 
morphological types. For example, iReddv fc Steidell {2009) 
found that a = —1.1 locally and increases up to ~ —1.7 at 
high redshift (from z=3 to 6). From our numerical exper- 
iments this case corresponds to P ~ 0.3. We can exclude 
all the cases with /3 > 0.4 and those with j3 < —0.2 since 
they predict values of a which are very much outside the ob- 
served range, and also a CSFR very much at variance with 
observations. 

In Figure [6] we show the case /? = 0.4 (top panel) 
which is similar to the suggested observational variation 
l|Reddv fc Steidell 120091 ) . and the case /3 = -0.2 (lower 
panel). Therefore, although we do not exclude a variation 
of the slope of the LFs of galaxies, we can conclude that the 
cases of a variations are providing a fit to the CSFR which is 
worse than the case in which there is number density varia- 
tion, as shown in Figure [4] Finally, we run several models by 
varying both <j>* and a and we did not find any combination 
of parameters which can fit the observed CSFR better than 
our best case with <j>* only variation. In Figure [7| we show 
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Figure 7. Comparison between our best model with NDE (blue 
solid line) and a model with both NDE (the same as the best 
model) and a evolution (green dashed line). Data points are the 
same described in Figure [3] 



a comparison between a model with the same variation of 
<f>* as our best model (see section 13 . 1 . 2 f) and a model with 
the same (f>* variation plus a variation of a, obtained with 
/3 = 0.3 for ellipticals and /3 = 0.2 for spirals and irregu- 
lars. As one can see, the agreement between the model with 
variable a and 4>* does not fit the data as well as our best 
model. 



3.1.4 Comparison with other determinations 

In Figure [8] we show the results of several determinations 
of the CSFR that we compare to our model predictions, in 
particular with the PLE and with our best model with NDE. 
The determinations of the CSFR we use for this comparison 



• \Steidel et all \l999\) : 

It is based on a sample of star-forming galaxies at 2 > 4. 
The correspondin g CSFR based on th e ir dat a has been then 
parametrized by iPorciani fc Madaul (|200ll ) (their model 
SF2L 

• I Cote et aij \200A) : 

Combining the data of the Two Micron All Sky Survey 
(2MASS) Extended Source Catalog and the 2dF Galaxy 
Redshift Survey, they produce an infrared selected galaxy 
catalogue. They use it to estimate the galaxy LF and to 
infer the total mass fraction in stars. Then they us e their 
results together with the data of ISteidel et alj l|l999l ) to ob- 
tain a parametric fit of the CSFR. Their param etric form 
has also been used b y Hopkin s fc Beacoml |2006). 

• \Porciani & Madaul \200j) : 



Figure 8. A compilation of determinations of the CSFH, com- 
pared to the results of this work. Mo dels are from: |S tcidcl et al. 
l|l999h magenta long da sh-dotted l i ne. [ Porciani fc Madaul l|200lh 
green short-dashed line. lCole et al. I ll200ll) dark gree n dotte d line, 
iMenci et all (|2004h blue solid line. IStrolger et al cyan 
dot short-dashed line, this work (PLE) black short-dashed long- 
dashed line, this work (NDE) turquoise long dashed line. Data 
point are the same as Figure [3] 



the hypothesis that the rate of GRBs traces the global star 
form ation history of the Universe. 
• \Strolser et al. 

It is a model bas ed on a modified ver sion of the parametric 
form suggested bv lMadau et al.l (|l998l ). taking into account 
dust extinction. Model parameters are deter mined fitting the 
collect ion of measurements of the CSFR of iGiavalisco et all 
(2004). 



Menci et ail \2004f) : 



Here we refer to their model SF3. It is a parameterization 
of the CSFR based on the data collected by the Burst and 
Transient Source Experiment (BATSE). It is adopted to test 



Prediction from a semi-analytic model for galaxy forma- 
tion in the framework of the hierarchical clustering scenario. 

For further information we address the reader to the 
above mentioned papers. 

From Figure [8] it is possible to see that our model 
with number density evolution shows a complete agree- 
ment with the semi analytical model of IMenci et all (2004) 
and, in general, shows a better agree ment with th e data , 
if c ompared to the determin ations of ISteidel et al.l (Il999l ) 
andU orciam fc Madaul l|200ll) . For what concerns the PLE 
model, this presents lower values than the other CSFRs, in 
the redshift range [0 — 6]. On the other hand, for z > 6 it 
exceeds all the other predictions considered here, as we have 
discussed already. 



3.1.5 The cosmic ISM mean metallicity 

Here we compute the evolution of the cosmic, luminosity 
weighted, mean interstellar metallicity of the galaxies of 
different morphological type. With this definition, we in- 
dicate the mean metallicity of the gas from which stars are 
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Figure 9. The cosmic luminosity weighted mean metallicity rela- 
tive to the Sun in the case of pure luminosity evolution. Magenta 
solid line: total; red dotted line: spirals; green long-dashed line: 
ellipticals; blue dashed line: irregulars. All the val ues are normal- 
ized t o the solar metallicty Zq = 0.0134, from Asp lund et al] 
j200Sh . 
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Figure 10. The cosmic luminosity weighted mean metallicity rel- 
ative to the Sun in the case of number density evolution. Magenta 
solid line: total; red dotted line: spirals; green long dashed line: 
ellipticals; blue dashed line: irregulars. All the val ues are normal- 
ized t o the solar metallicty Zq = 0.0134, from Asplu nd et al] 
J2009h . 



born, in a unita r y vo lume of the Universe. As suggested by 
iKulkarni fc Fall (2002), we compute it through the following 
expression: 



/ Zi(Lj)Li<f>i(Li) dLj 
J2t J Li^i(Li) dLi 



(15) 



where Zi(L) is the average interstellar metallicity in a galaxy 
of luminosity Li , at any given cosmic time, and of the i-th 
morphological type and $i(Li) is the luminosity function of 
the i-th morphological type. 

We compute this quantity both in the case of PLE and 
NDE, using the metallicity evolutions obtained with the 
chemical evolution models, together with the local B-band 
luminosities obtained with the spectrophotometric code. 
The parameters of the local luminosity functions, are from 
Mar zke et al.l (1 19981 ). as in the previous paragraphs. 

From Figure [9] it is clear that in the PLE scenario el- 
lipticals dominate the mean galactic metallicity throughout 
the whole cosmic time. This is due to their higher average 
metallicity compared to spirals and irregulars and to the 
constant number density of all galaxies. The spiral galax- 
ies are the second most important contributors to the cos- 
mic chemical enrichment, whereas irregulars give a negligi- 
ble c ontribution. We can see tha t our model is in agreement 
with ICalura fc Matteuccil (|2004l ) who found that for z > 6 
the main metal producers are ellipticals whereas for z < 2 
the main sources of metals production are spirals. It is worth 
noting that, the decrease of the luminosity weighted metal- 
licity of ellipticals is due to the fact that, at z ~ 6 in our 
model, ellipticals stop forming stars and from that moment 
on their B-band luminosity abruptly decreases. Also spirals 
decrease their B-band luminosity but mildly since their star 



formation lasts for the whole cosmic time. We note that in 
this case two distinct effects are involved: the metallicity 
which is increasing and the luminosity which is decreasing 
with decreasing z. 

In Figure [10] we show the same as in Figure [9] but under 
the hypothesis of number density evolution. In particular, 
we show the results of our best model, as described in sec- 
tion [3T]2] The most evident fact in this case is that spirals, 
instead of ellipticals, are the main contributors to the to- 
tal cosmic luminosity weighted mean metallicity. However, 
we need further investigation to better interpret this result. 
It has to be considered that this quantity cannot be com- 
pared directly with observations, since it is a cosmic aver- 
age, while measurements are performed in single objects. 
Another important thing we note is that, in this case the 
cosmic luminosity weighted mean metallicity of ellipticals 
is increasing throughout the whole cosmic time. This can 
be explained by the fact that the decreasing of the B-band 
luminosity is compensated by the increase of their number 
density at lower redshift. We also find that our present time 
value of the cosmic luminosity weighted mean metallicity of 
the ISM in all galaxies and in each galaxy formation scenario 
is ~ 2.4 Zq, a value that is sl i ghtly higher than 1.3 Zq as 
found bv lCalura fc Matteuccil (|2004h . This i s partly due to 
the fa ct that we adopt the solar metallicity of lAsplund et al] 
(2009), that is lower than the Zq adopted in the previous 
paper. 

This average metallicity is not the global average metal- 
licity of the Universe since we did not consider the metallic- 
ity of the intracluster and intergalactic medium. 
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4 DISCUSSION AND CONCLUSIONS 

In this paper we computed the CSFR by means of de- 
tailed chemical and spectro-photometric models of galaxies 
of different morphological type (ellipticals, spirals, irregu- 
lars). These models well reproduce the main properties of 
local galaxies and have been tes ted in the past o n a variety 
of observational constraints (see lMatteu cci 2001). Models of 
galactic chemical evolution of galaxies contain a few free pa- 
rameters such as the efficiency of star formation, the slope of 
the IMF in the case of external galaxies, the fraction of bi- 
nary systems in the IMF giving rise to SNe la, the timescale 
for gas accretion and the efficiency of possible galactic winds. 
All these parameters should be fixed by reproducing the 
largest as possible number of observational constraints, so 
that the number of well fitted features outnumbers the num- 
ber of free parameters. For the Milky Way we have many ob- 
servational features to reproduce, such as the G-dwarf metal- 
licity distribution, the solar abundances, the present time 
amount of gas, the present time SN rates, the star formation 
and gas infall rates, the abundance patterns in stars ([X/Fe] 
vs. [Fe/H], where X represents the abundance of a chem- 
ical elements). Then one should reproduce the abundance 
gradients of chemical abundances, star formation rate and 
gas distribution along the disk. For external galaxies, such 
as ellipticals, the number of constraints is less than for the 
Milky Way but still they are enough to counterbalance the 
number of free parameters involved in the models. The ob- 
servational features to reproduce are: the Color-Magnitude 
diagram, the average metallicity (measured by Mg and Fe) 
of the stars, the abundance gradients in the stellar popula- 
tions as a function of the galactocentric distance, the high 
[a/Fe] ratios measured in the dominant stellar population 
in the central part of these galaxies, the mass-metalicity re- 
lation and the [a/Fe] vs. mass relation, the Type la SN 
rate at the present time. This last features implies that 
the more massive galaxies show higher central [a/Fe] ra- 
tios than less massive ones. This is explained at the moment 
only by chemical evolution models of the knid presented 
here. Classical hiera rchical models predict the contrary (see 
iThomas fc al.ll2002l'). while mo r e recent semi-analyt ical mod- 
els (e.g. iDe Lucia et all 120061 ; iPipino et all 120091 ) improve 
the situation but they still do not reproduce entirely this 
relation. In order to obtain an increase of the [a/Fe] ratio 
with galactic mass is enough to assume an increasing star 
formation efficiency with galactic mass. However, we cannot 
exclude a priori other solutions, such as a variable IMF from 
galaxy to gala xy, but this solution creates other problems 
(see iMatteucci et al.lll99"H ). Also the models for irregulars 
have been tested and they reproduce the metallicity, the 
amount of gas, the [C,N,0/Fe] vs. [Fe/H] relations, the SN 
rates for an average irregular galaxy. Therefore, we felt rea- 
sonably confident that the star formation histories suggested 
by these galaxy models could be used to compute the CSFR. 
We then adopted the histories of star formation predicted 
by such models and convolved them with the luminosity 
functions of galaxies. These LFs have been constrained to 
reproduce the observed local LFs. We first computed the 
CSFR under the assumption of pure luminosity evolution of 
galaxies, namely we did assume that the main parameters of 
the LF do not vary with redshift, in particular the number 
density of galaxies and the slope of the LF. By doing that, 



we ob tained results similar to those of ICalura fc Matteuccil 
(2003). Then we computed the CSFR by assuming number 
density and slope evolution. 

Our results can be summarized as follows: 

• In the case of pure luminosity evolution of galaxies, the 
CSFR traces the star formation histories of galaxies. In the 
framework of our galaxy evolution models, ellipticals form 
at very high redshift in a very short interval of time and 
then evolve passively for the remaining galactic age. The 
spirals and even more the irregulars, instead, are forming 
stars continuously and at a much lower rate for the whole 
galactic lifetime. In this situation, we predict that the ellip- 
ticals are the dominant source of the cosmic star formation 
at very high redshift (from z — 6 to z — 10), whereas spi- 
rals start to dominate for z < 6. On the other hand the 
contribution of irregulars to both the CSFR and the cosmic 
chemical enrichment is found to be negligible. By compar- 
ing these model results with the available data for CSFR, 
we conclude that our models overpredict the CSFR at high 
redshift and underpredict it at low redshift. However, most 
of models for CSFR in the literature also underpredict this 
quantity at low redshift. For the high redshift situation we 
should take into account that the data for z > 6 are still very 
uncertain, especially because of the uncertain dust correc- 
tions and assumed IMF. Therefore, a conclusive statement 
on this point would be premature. For the low redshift we 
probably underestimate the amount of cosmic star forma- 
tion, since we did not include massive starburst galaxies at 
intermediate and low redshift. 

• In the case of number density evolution, by assuming 
that ellipticals were less at high redshift and their num- 
ber increased at lower redshifts, and that the contrary hap- 
pened to spirals, we obtain a CSFR in much better agree- 
ment with observations. This result could imply that many 
ellipticals formed by merging of sp irals, as originally sug- 
gested bv lToomre fc Toomrel l|l972h . In particular, our pre- 
dicted number density evolution fo r ellipticals impli e s 4>* = 
0o (1 + z)~ 8 in agreement with iTotani fc Yoshiil (|l99St ). 
wher eas for spirals impl i es 4>* = (j>o (l + z) in ag r eemen t 
with Ide Lapparent et al.l (120041 1 and iBoissier et al.l l|2010f) . 
Possible criticism to the scenario of the dry mergers of spi- 
rals to form ellipticals is given by the fact that a dry merg- 
ing process seems unrealistic in the light of the difficulty of 
reproducing the observed local chemical properties of ellip- 
ticals. In fact, it is well known that the chemical properties 
of spirals and ellipticals are different (e.g. different average 
metallicities and [a/Fe] ratios). By comparing our predicted 
CSFR with other determinations of the CSFR we found that 
it is similar to the CSFR of the hierarchical model of Menci 
& al. (2004). 

• We tested also the case of a variation of the slope of 
the luminosity function with redshift and concluded that a 
variation of the slope similar to what is suggested by ob- 
servations is not excluded, although it does not provide an 
improvement relative to pure luminosity evolution case and 
it predicts a CSFR in worse agreement with the data than 
the case of only galaxy number density variation. Also cases 
with both NDE and slope evolution at the same time do not 
provide better fits. 

• We computed the cosmic, luminosity weighted, metal- 
licity in the gas of galaxies and its evolution with redshift. 
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We found that, in the case of pure luminosity evolution, the 
ellipticals dominate the cosmic metal enrichment over the 
whole lifetime of the Universe. On the other hand, in the 
case of number density evolution, are the spirals to domi- 
nate the cosmic metal enrichment of the Universe. This is 
potentially an important result but it needs to be confirmed 
by future studies. 
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